Hypothesis testing for an entangled state produced 
by spontaneous parametric down conversion 



O 
O 

(N 

bo: 

< 
(N 

(N 
> 

(N 

m 
O 

O 

Oh! 

H— > 

CT 1 - 



X 



Masahito Hayashi*d, Bao-Sen Shi*, Akihisa Tomita*'*, 
Keiji Matsumoto*'§, Yoshiyuki Tsuda^, and Yun-Kun Jiang* 
* ERATO Quantum Computation and Information Project, 
Japan Science and Technology Agency (JST), Tokyo 113-0033, Japan 
' Superrobust Computation Project, Information Science and Technology Strategic Core (21st Century COE by MEXT) 
Graduate School of Information Science and Technology, The University of Tokyo 
7-3-1, Hongo, Bunkyo-ku, Tokyo, 113-0033, Japan 
^Fundamental Research Laboratories, NEC, Tsukuba 305-8501, Japan 
^National Institute of Informatics, Hitotsubashi, Chiyoda-ku, Tokyo 101-8430, Japan 
1 COE, Chuo University, Kasuga, Bunkyo-ku, Tokyo 112-8551, Japan 

Generation and characterization of entanglement are crucial tasks in quantum information pro- 
cessing. A hypothesis testing scheme for entanglement has been formulated. Three designs were 
proposed to test the entangled photon states created by the spontaneous parametric down con- 
version. The time allocations between the measurement vectors were designed to consider the 
anisotropic deviation of the generated photon states from the maximally entangled states. The 
designs were evaluated in terms of the p-value based on the observed data. It has been experimen- 
tally demonstrated that the optimal time allocation between the coincidence and anti-coincidence 
measurement vectors improves the entanglement test. A further improvement is also experimentally 
demonstrated by optimizing the time allocation between the anti-coincidence vectors. Analysis on 
the data obtained in the experiment verified the advantage of the entanglement test designed by the 
optimal time allocation. 

PACS numbers: 03.65.Wj,42.50.-p,03.65.Ud 



I. INTRODUCTION 

The concept of entanglement has been thought to be 
the heart of quantum mechanics. The seminal experi- 
ment by Aspect et al. PJ has proved the 'spooky' non- 
local action of quantum mechanics by observing violation 
of Bell inequality with entangled photon pairs. Re- 
cently, entanglement has been also recognized as an im- 
portant resource for quantum information processing, ex- 
plicitly or implicitly. For example, entanglement provides 
an exponential speed-up in some computational tasks Q , 
and unconditional security in cryptographic communica- 
tions 0. A hidden entanglement between the legitimate 
parties guarantees the security in BB84 quantum cryp- 
tographic protocol 0, Q . Quantum communication be- 
tween arbitrarily distant parties has been shown to be 
possible by a quantum repeater [f| based on quantum 
teleportation pfl • Practical realization of entangled states 
is therefore one of the most important issues in the quan- 
tum information technology. 

The practical implementation raises a problem to ver- 
ify the amount of entanglement. A quantum informa- 
tion protocol requires a minimum entanglement. It is, 
however, not always satisfied in actual experiments. Un- 
avoidable imperfections will limit the entanglement in 
generation process. Moreover, decoherence and dissipa- 
tion due to the coupling with the environment will de- 
grade the entanglement during the processing. Therefore, 
it is a crucial issue to characterize the entanglement of 
the generated (or stored) states to guarantee the success- 
ful quantum information processing. For this purpose, 



quantum state estimation and Quantum state tomogra- 
phy are known as a method of identifying the unknown 
state d, U ^3 ■ Quantum state tomography ^ nas re ~ 
cently applied to obtain full information of the 4x4 two- 
particle density matrix from the coincidence counts of 16 
combinations of measurement |12(. However, characteri- 
zation is not the goal of an experiment, but only a part 
of preparation. It is thus favorable to reduce the time 
for characterization and the number of consumed parti- 
cles as possible. An entanglement test, which should be 
simpler than the full characterization, works well in most 
applications, because we only need to know whether the 
states are sufficiently entangled or not. We can reduce 
the resources for characterization with the entanglement 
test. Barbieri et al. 01 introduced an entanglement 
witness to test the entanglement of polarized entangled 
photon pairs. Tsuda et al. |l4j studied the optimization 
problem on the entanglement tests with the mathemat- 
ical statistics in the POVM framework. Hayashi et al. 
[l5j treated this optimization problem in the framework 
of Poisson distribution, which describes the stochastic be- 
havior of the measurement outcomes on the two-photon 
pairs generated by spontaneous parametric down conver- 
sion (SPDC). 

In this paper, we will apply the experimental designs 
proposed in ref. |15| to test the polarization entangled 
two-photon pairs generated by SPDC. The two-photon 
states can be characterized by the correlation of photon 
detection events in several measurement bases. In exper- 
iments, the correlation is measured by coincidence counts 
of photon detections on selected polarizations. The co- 
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incidence counts on a combination of the two from hori- 
zontal (H), vertical (V), 45° linear (X), 135° linear (£)), 
clock- wise circular (i?), and anti-clock- wise circular (L) 
polarizations defines a measurement vector. We here use 
the same representations for the measurement vectors as 
the state vectors, such as \HH). If the state is close 
to |$(+)> = -X={\HH) + \VV)), the coincidence counts 
on the vectors \HH), \VV), \DD), \XX), \RL), and \LR) 
yield the maximum values, whereas the counts on the vec- 
tors \HV},\VH),\DX),\XD) i \RR), and \LL) take the 
minimum values. We will refer to the former vectors as 
the coincidence vectors, and to the latter as the anti- 
coincidence vectors. The ratio of the minimum counts 
to the maximum counts measures the degree of entan- 
glement. In the following sections, we formulate the hy- 
pothesis testing of entanglement in the view of statistics. 
We then improve the test by optimizing the allocation 
on the measurement time for each measurement, consid- 
ering that the counts on the anti-coincidence vectors are 
much smaller than those on the coincidence vectors. 



The test can be further improved, if we utilize the 
knowledge on the tendency of the entanglement degrada- 
tion. In general, the error from the maximally entangled 
states can be anisotropic, which reflects the generation 
process of the states. We can improve the sensitivity to 
the entanglement degradation by focusing the measure- 
ment on the expected error directions. In the present 
experiment, we generated polarization entangled photon 
pairs from a stack of two type-I phase matched nonlin- 
ear crystals [l6lll7j. where one nonlinear crystal gener- 
ates a photon pair polarized in the horizontal direction 
(\HH)), and the other generates a photon pair polarized 
in the vertical direction If the two pairs are in- 

distinguishable, the generated photons are entangled in 
a two-photon state -^(\HH) + cxp[iS]\VV)) . Otherwise, 
the state will be a mixture of HH pairs and VV pairs. 
The quantum state tomography has shown that the only 
HHHH, WW, VVHH, HHVV elements are domi- 
nant p"?||. which implies that the density matrix can be 
approximated by a classical mixture of |$( + )) ($( + ) | and 
|$(-))($(-)|. We can improve the entanglement test on 
the basis of this property, as described in the following 
sections. 



The construction of this article is following. Section ITT1 
gives the mathematical formulation concerning statisti- 
cal hypothesis testing. Section ITTT1 defines the hypothesis 
scheme for the entanglement of the two-photon states 
generated by SPDC. Sections llVl - IVIll describe testing 
methods as well as design of experiment. Section lYlIIl ex- 
amines the experimental aspects of the hypothesis testing 
on the entanglement. The designs on the time allocation 
are evaluated by the experimental data. 



II. HYPOTHESIS TESTING FOR 
PROBABILITY DISTRIBUTIONS 

A. Formulation 

In this section, we review the fundamental knowledge 
of hypothesis testing for probability distributions [lg. 
Suppose that a random variable X is distributed accord- 
ing to a probability measure Pg identified by an unknown 
parameter 9, and that the unknown parameter 9 belongs 
to one of the mutually disjoint sets ©o and 0i. When the 
task is to guarantee that the true parameter 9 belongs to 
the set 0i with a certain significance, we choose the null 
hypothesis Hq and the alternative hypothesis Hi as 



Hq : 9 £ ©o vs. Hi : 9 S X . 



(1) 



The task is then described by a test, where we decide to 
accept the hypothesis Hi by rejecting the null hypothesis 
Hq with confidence. We make no decision when we can 
not reject the null hypothesis. The test based is charac- 
terized as a function 4>{x) taking a value in {0, 1}; we can 
reject Hq if 4>{x) = 1, but can not reject it if 4>{x) = 0. 
The test can be also described by the rejection region 
defined by {x\<fi(x) = 1}, whether the data x falls into 
the rejection region or not. The function 4>(x) should be 
designed properly according to the problem to provide 
an appropriate test. 

For example, suppose that mutually independent data 
xi,--- ,x n obey a normal distribution with unknown 
mean /i and known variance erg 



: exp 



We introduce a hypothesis testing 



(2) 



H : n E < no} vs. Hi : fie {fj,\fi > fio} (3) 

to decide the mean is larger than a value /iq with a confi- 
dence level a. Then the function cf> of the average on the 
data is defined by x = (1/n) J2™ = 



fa) 



1, if ^/n(x - no)/<Jo > z a 
0, if y/n(x - Ho)/<Jo < z a , 



where z a is given by 
$(z) 



1 



2tt 



e 2 dy. 



The rejection region is provided as 

{x\y/n(x - ^ )/(70 > Zq] 

in this test. 



(4) 

(5) 
(6) 

(7) 
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B. p- values 

In a usual hypothesis testing, we fix our test before ap- 
plying it to the data. Sometimes, however, we compare 
tests in a class T by the minimum risk probability to re- 
ject the hypothesis Hq with given data. This probability 
is called the p- value, which depends on the observed data 
x as well as the subset Go to be rejected. For a given class 
T of tests, the p- value is defined as 



pval := min max Pg (</>). 



(8) 



In the hypothesis test for a normal distribution defined 
previously (Eqs. ©-©), the p- value is given by 



pval = <&(y/n(x - fx )/a ). 



(9) 



The concept of p- value is useful for comparison of several 
designs of experiment. 



C. Likelihood Test 

In mathematical statistics, the likelihood ratio tests is 
often used as a class of standard tests^8|- Likelihood is 
a function of unknown parameter 9 defined as the prob- 
ability distribution Pg{x) for given data x. When both 

00 and 9i consist of single elements as Qq = {9q} and 

01 = the rejection region of the likelihood ratio 
test is defined as 



Pb (x) 



9iO) 



< r 



where r is a constant, and the ratio Pg (x)/ Pg 1 (x) is 
called the likelihood ratio. In the general case, i.e., the 
cases where ©o or <di has at least two elements, the like- 
lihood ratio test is given by its rejection region: 



su Peee p e( x ) 



< r 



In the hypothesis test for a normal distribution defined 
previously (Eqs. Q-®); logarithm of the likelihood ra- 
tio is given by 



log 



Pp,=s{x) 2a, 



- mo) 2 - Yl 



n(x - fi ) 2 



(xi - ^) 2 } 
(10) 



Therefore, the rejection region of the likelihood ratio test 
coincide with Eq. 0, if we set r = exp[— z^/2]. 



pairs generated by SPDC. The two-photon state is de- 
scribed by a density matrix a. We assume two-photon 
generation process to be identical but individual. Here 
we measure the entanglement by the fidelity between the 
generated state a and the maximally entangled target 
state |$( + )): 



(11) 



The purpose of the test is to guarantee, with a certain 
significance, that the state is sufficiently close to the max- 
imally entangled state. For this purpose, the hypothesis 
that the fidelity F is less than a threshold Fq should be 
disproved with a small error probability. In mathematical 
statistics, the above situation is formulated as hypothesis 
testing; we introduce the null hypothesis Hq that entan- 
glement is not enough and the alternative Hi that the 
entanglement is enough: 



H ■ F < F vs. H x : F > F , 



(12) 



with a threshold Fo. 

The coincidence counts of photon pairs generated in 
the SPDC experiments can be assumed to be a random 
variable independently distributed according to a Poisson 
distribution. In the following, a symbol labeled by a pair 
(x, y) represents a random variable or parameter related 
to the measurement vectors \x,y). When the dark count 
is negligible, the number of detection events (i.e., coinci- 
dence counts) n xy on the vectors \x,y) is a random vari- 
able according to the Poisson distribution Poi(\fi xy t xy ) 
of mean \p, xv t xyi where 

• A is a known constant related to the photon detec- 
tion rate, determined from the averaged photon- 
pair-generation rate and the detection efficiency, 

• l^xy = (x, y|c|a;, y) is an unknown constant, 

• t xy is a known constant of the time for detection. 
The probability function of n xy is 



n xy = Poi(X/j, xy t xy ) = e^p{-X^i xy t xy ) 



(X[i X yt xy ^ Xy 



Tl xy ! 



Because the detections at different times are mutually in- 
dependent, n xy is independent of n x i y : {x ^ x' or y ^ y 1 ) 
In this paper, we discuss the quantum hypothesis test in 
under the above assumptions, whereas Usami et al 
discussed the state estimation under the same assump- 
tions. 



III. HYPOTHESIS TESTING SCHEME FOR 
ENTANGLEMENT IN SPDC EXPERIMENTS 



B. Modified visibility 



A. Formalization 

This section introduces the hypothesis test for entan- 
glement. We consider the entanglement of two-photon 



Visibility of the two-photon interference is an indica- 
tor of entanglement commonly used in the experiments. 
The two-photon interference fringe is obtained by the 
measurement of coincidence counts on the vector \x,y), 
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where the vector \y) is rotated along a great circle on 
the Poincare sphere with a fixed vector |x). The visibil- 
ity is calculated with the maximum and minimum num- 
ber of coincidence counts, n max and n m i n , as the ratio 
{n max - n min )/(n max + n min ). We need to make the 
measurement with at least two fixed vectors \x) in order 
to exclude the possibility of the classical correlation. We 
may choose the two vectors \H) and \D) as \x), for ex- 
ample. However, our decision will contain a bias, if we 
measure the coincidence counts only with two fixed vec- 
tors. The bias in the measurement emerges as the fact 
that the visibility reflects not only the fidelity but also 
the direction of the deviation of the given state from the 
maximally entangled target state. Hence, we cannot es- 
timate the fidelity in a statistically proper way from the 
visibility. 

In order to remove the bias based on such a di- 
rection, we propose to measure the counts on the 
coincidence vectors \HH), \VV), \DD), \XX), \RL), and 
| Li?), and the counts of the anti-coincidence vectors 
\HV),\VH),\DX),\XD),\RR), and \LL). The former 
corresponds to the maximum coincidence counts in the 
two-photon interference, and the latter does to the min- 
imum. Since the equations 

\HH)(HH\ + \VV)(VV\ + \DD)(DD\ 
+ \XX)(XX\ + \RL)(RL\ + \LR)(LR\ 
=2|$(+))($(+)| +/ (13) 

and 

\HV)(HV\ + \VH)(VH\ + \XD)(XD\ 
+ \DX)(DX\ + \RR)(RR\ + \LL)(LL\ 
=2/-2|$( +) )($(+)| (14) 

hold, In this paper, we call this proposed method the 
modified visibility method. Using this method, we can 
test the fidelity between the maximally entangled state 
|$(+)) ($(+) | and the given state a, using the total num- 
ber of counts of the coincidence events (the total count 
on coincidence event) n\ and the total number of counts 
of the anti-coincidence events (the total count on anti- 
coincidence events) ri2 obtained by measuring on all the 
vectors with the time 4r. holds, we can estimate the 
fidelity by measuring the sum of the counts of the follow- 
ing vectors: \HH)AVV\\DD),\XX),\RL), and \LR), 
when A is known [la. fl4| . This is because the sum 
ni := n H H + n V v + n DD + n X x + n RL + n LR obeys the 
Poisson distribution with the expectation value A 1+ 6 t\ , 
where the measurement time for each vector is We 
call these vectors the coincidence vectors because these 
correspond to the coincidence events. 

However, since the parameter A is usually unknown, we 
need to perform another measurement on different vec- 
tors to obtain additional information, also holds, we can 
estimate the fidelity by measuring the sum of the counts 
of the following vectors: \HV), \VH), \DX), \XD), \RR), 
and ILL). The sum n 2 := nnv + n VH + nux + nxD + 



nRR + riLL obeys the Poisson distribution Poi(A 6 t 2 ), 
where the measurement time for each vector is ^ . Com- 
bining the two measurements, we can estimate the fi- 
delity without the knowledge of A. We call these vectors 
the anti-coincidence vectors because these correspond to 
the anti-coincidence events. 

We can also consider different type of measurement on 
A. If we prepare our device to detect all photons, the 
detected number 713 obeys the distribution Poi(At 3 ) with 
the measurement time £3. We will refer to it as the total 
flux measurement. In the following, we consider the best 
time allocation for estimation and test on the fidelity, 
by applying methods of mathematical statistics. We will 
assume that A is known or estimated from the detected 
number 77,3. 



IV. MODIFICATION OF VISIBILITY 

The total count on coincidence events n\ obeys 
Poi(A^rtii), and the count on total anti-coincidence 
events 712 obeys the distribution Poi(A 2 ~ 2F £). These ex- 
pectation values fi\ and ^2 are given as (i± = A 2 ^ 1 ?; 



and \ii 



A 



2-2F 
12 



t. Since the ratio 



Ml+A»2 



is equal to 



|(1 — F), we can estimate the fidelity using the ratio 



" 2 as F{n\,n-2) 
totically equal to 



712 



2 ni+n 2 



Its variance is asymp- 



A( 



3(2F+1) ' 3(2-2F) < 



(2F+1)(2-2F) 
At ' 



(15) 



Hence, similarly to the visibility, we can check the fidelity 
by using this ratio. 

Indeed, when we consider the distribution under the 
condition that the total count n\ + ni is fixed to n, 
the random variable obeys the binomial distribu- 
tion with the average value |(1 — F)n. Hence, we 
can apply the likelihood test of the binomial distribu- 
tion. In this case, by the approximation to the nor- 
mal distribution, the likelihood test with the risk prob- 
ability a is almost equal to the test with the rejection 
region concerning the null hypothesis Ho : F < Fq: 

{(ni,n 2 )\-%- < | ( 1 - F ) + S-W ^&r* }. 



where $(a) := J" 
of tests is $( 



m-Fo) + $-\a) Jl 9{ni+n2) 

e~~ dx. The p- value of this kind 



1 

/27T ' 



K 2 (2F + l)-n 1 (2-2F ) 



V(«i+«2)(2Fo + l)(2-2F ) 



)• 



V. DESIGN I (A: UNKNOWN, ONE-STAGE) 

In this section, we consider the problem of test- 
ing the fidelity between the maximally entangled state 
|$(+))($(+) I an d the given state a by performing three 
kinds of measurement, coincidence, anti-coincidence, and 
total flux, with the times t\,t% and <3, respectively. The 
data (ni, ri2, 713) obeys the multinomial Poisson distri- 
bution Poi(A 2F r +1 ti, A 2 ~ 2F 7;2, Ar;3) with the assumption 
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that the parameter A is unknown. In this problem, it is 
natural to assume that we can select the time allocation 
with the constraint for the total time U + t 2 + t^ = t. 

As is shown in Hayashi et al when F\ := 

0.899519 < F < 1, the minimum variance of the estima- 
tor is asymptotically equal to ( 1 ~- F )( v/ |+ v/1 F ^ , which 
is attained by the optimal time allocation t\ = 0, 



s/S 



-t, ta 



=i. Otherwise, it 



(2F+l)(l-F)(V2F+l + \/2-2F) 



(V3+VT=F) ' 

is asymptotically equal to 

which is attained by the optimal time allocation t\ = 
ta = ,,.^^+1 wmt ta = 0. This op- 



3At 



V2F+1 + V2-2F"' "~ z \/2F+l + V2-2F*' ^ 3 

timal asymptotic variance is much better than that ob- 
tained by the modified visibility method. 

Since the optimal allocation depends only on the true 
parameter F, it is suitable to choose the optimal time 
allocation at the threshold Fq for testing whether the 
fidelity is greater than the threshold Fq. In the both 
cases, when we perform the optimal allocation, the data 
obeys the binomial Poisson distribution. Hence, similarly 
to the modified visibility, we can check the fidelity by 
applying the likelihood test of binomial distribution to 
the ratio between two kinds of counts. 

In the former case, the likelihood test with the risk 
probability a is almost equal to the test with the rejec- 
tion region Hq : F < F : {(n 2 ,n 3 )|— ^ < / y i ~ F ° + 

}. The p- value of this kind of tests 
). In the later case, the likelihood 



y/3 + y/l-F V 112+713 



V("2+™3)v /ITr Po'v / 3 / 

test with the risk probability a is almost equal to the test 
with the rejection region concerning the null hypothe- 
sis H : F < F : {(n l ,n 2 )\-, 

* _1 (c0 / V2-2FqV^F^H 



1 n 2 +ni 



< 



V2Fo + l + V2-2F 



V2Fo + l + x/2-2F 

kind of tests is $( 



711+T12 



ri2 \Z2F Q + l-7ti V2-2F 



V(™i+™2)\ / 2F + l % /2 rr 2^ 



The p-value of this 
)■ 



VI. DESIGN II (A: KNOWN, ONE-STAGE) 

In this section, we consider the case where A is known. 
In this case, when F > j, the optimal time allocation is 
t 2 = t, t\ = i;3 = 0. That is, the count on anti-coincidence 
(ti = 0; t 2 = f) is better than the count on coincidence 
(<i = t]t% = 0). In fact, Barbieri et aL[13j measured 
the sum of the counts on the anti-coincidence vectors 
\HV),\VH),\DX),\XD),\RR),\LL) to realize the en- 
tanglement witness in their experiment. In this case, the 
count on anti-coincidence n 2 obeys the Poisson distribu- 
tion with the average X^-^-t, the fidelity is estimated by 

1 — 3 4| L . Its variance is 3 Z. ■ Then, the likelihood 
test with the risk probability a of the Poisson distribu- 
tion is almost equal to the test with the rejection region: 

{n 2 \n 2 < A(1 3 fo) t + $- 1 (a)y^^2l*} concerning the 
null hypothesis Hq : F < Fq. The p-value of likelihood 
tests is $(" 2 7 A(1 - Fo)t/3 ). 



When F < j, the optimal time allocation is t\ = t, 
t 2 = h = 0. The fidelity is estimated by ^ - \. Its 

variance is ^^x^ . Then, the likelihood test with the 
risk probability a of the Poisson distribution is almost 
equal to the test with the rejection region: {ni|ni > 

A i ±f^i+$- 1 (l-a)y / A 1± f Fa t} concerning the null hy- 
pothesis Hq : F < Fq. The p-value of likelihood tests is 

^/■ -m+A(l+2F )t/6 N 
1 y/X(l+2F )t/6 >' 

VII. DESIGN III (A: KNOWN, TWO-STAGE) 

Next, for a further improvement, we minimize the 
variance by optimizing the time allocation tnv, tvH, 
tDX) txD, tRR, and t^j, between the anti-coincidence 
vectors B = \\HV), \VH), \DX), \XD), \RR), \LL)}, 
respectively, under the restriction of the total measure- 
ment time: t HV + t VH + t DX + txD + tRR + t LL = t. 
The number of the counts n xy obeys Poisson distribu- 
tion Poi(Xfj, xy t xy ) with the unknown parameter /i xy — 
(xaHbI^IxaUb) ■ the minimum value of estimation of the 
fidelity F = 1 - | T,( x , y )eB tow is 

(E 

AXt 

which is attained by the optimal time allocation 



t x y — 



Yl<(x,y)eB y/Vx,y 



(16) 



(17) 



which is called Neyman allocation and is used in sampling 
design 19]. 

However, this optimal time allocation is not applicable 
in the experiment, because it depends on the unknown 
parameters [Iflv ■, MVff, HdX, Vxd, Hrr, and ull- In or- 
der to resolve this problem, we use a two-stage method, 
where the total measurement time t is divided into tf 
for the first stage and t s for the second stage under the 
condition of t = tf + t s . In the first stage, we mea- 
sure the counts on the anti-coincidence vectors for f//6 
and estimate the expectation value for Neyman alloca- 
tion on measurement time t s . In the second stage, we 
measure the counts on anti-coincidence vectors \xaVb) 
according to the estimated Neyman allocation. The two- 
stage method is formulated as follows. 

(i) The measurement time for each vector in the first 
stage is given by t//6 

(ii) In the second stage, we measure the counts on a vec- 
tor \xaVb) the measurement time t xy defined as 



txy 



xy 



J2(x,y)eB \A 



:(*"*/) 



' xy 



where m xy is the observed count in the first stage, 
(hi) Define £i xy and F as 



H>xy 



L xy 



At, 



F=l- 1 - 
2 



x.y)eB 



6 



where n x , y is the number of the counts on \xaVb) for i xy . 
Then, we can estimate the fidelity by F. 

We can test our hypothesis by using the likelihood test, 
which has the rejection region: 



sup Poi(/2)(n) < r sup Poi(/t)(n) 



where Wi := ^tt" an d cq := 1 — Fq. However, it is very 
difficult to choose r such that the likelihood test has a 
given risk probability a. In order to resolve this problem, 
Hayashi et al. 15] proposed a better method. For treat- 
ing this method, we number the anti-coincidence vectors 
by i = 1, . . . , 6. Then, the rejection region with the risk 
probability a is given by Br , where the set Br is defined 
by 



B R := 



and fii(R) are defined as follows: 



Hi(R) + (j,i{R) log 

w l c 



< 1 



(18) 



R ifi?<i? , 4 (19) 



^+f ll {R)\og WM Wi =R HR>R 0t (20) 
w M w M 



where w m 



log 



maxi Wi and i?o,i := 
The value R a is defined by 



Co 



min Zij(Ra) — <f> 1 (a), 



(21) 



where Zjj (i?) is given by 



Xj{R) 



if 



2x J {R)y i (R) 



> 1 



x i (R)y i (R)+x i (R)y j (R) 

■ r 2x l (R)y J (R) > i 

:Cj (i?)y J (i?)+ 2;j (i?)y I (i?) ^ 



(22) 



Otherwise, 



2(x i (R)x j (R)(y i (R) + yj (R)) - Xi (R) 2 yj (R) - Xj (R) 2 yi (R)) 



y/(xi(R) - x 3 (R))( yi (R) - y~<M^x~{Ejy~{Rf + Xj (R) yi (R)* - y^y^R^R) + Xj (R)) ' 



(23) 



r 



where x t (R) := jrfjfa 



land m {R) :=^p 



Concerning these tests, the p-value is calculated to 

$(-mmz^(i?'(n))), (24) 

where i?'(n) is defined as 

6 



V — ^ — = i. 



(25) 



VIII. ANALYSIS OF EXPERIMENTAL DATA 

The experimental set-up for the hypothesis testing is 
shown in Fig. 2] The nonlinear crystals (BBO), the op- 
tical axis of which were set to orthogonal to on another, 
were pumped by a pulsed UV light polarized in 45° di- 
rection to the optical axis of the crystals. One nonlinear 
crystal generates two photons polarized in the horizon- 
tal direction (\HH)) from the vertical component of the 
pump light, and the other generates ones polarized in the 



vertical direction (|VV)) from the horizontal component 
of the pump. The second harmonic of the mode-locked 
Ti:S laser light of about 100 fs duration and 150 mW 
average power was used to pump the nonlinear crystal. 
The wavelength of SPDC photons was thus 800 nm. The 
group velocity dispersion and birefringence in the crys- 
tal may differ the space-time position of the generated 
photons and make the two processes to be distinguished 
|17| . Fortunately, this timing information can be erased 
by compensation; the horizontal component of the pump 
pulse should arrive at the nonlinear crystals earlier than 
the vertical component. The compensation can be done 
by putting a set of birefringence plates (quartz) and a 
variable wave-plate before the crystals. We could con- 
trol the two photon state from highly entangled states to 
separable states by shifting the compensation from the 
optimal setting. 

The count on the vector \xaTJb) was measured by ad- 
justing the half wave plates (HWPs) and the quarter 
wave plates (QWPs) in Fig. ^ We accumulated the 
counts for one second, and recorded the counts every one 
second. Therefore, the time allocation of the measure- 
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FIG. 1: Schematic of the entangled photon pair generation 
by spontaneous parametric down conversion. Cascade of the 
nonlinear crystals (NLC) generate the photon pairs. Group 
velocity dispersion and birefringence in the NLCs are pre- 
compensated with quartz plates and a Bereck compensator. 
Two-photon states are analyzed with half wave plates (HWP) , 
quarter wave plates (QWP), and polarization beam splitters 
(PBS). Interference filters (IF) are placed before the single 
photon counting modules (SPCM). 



FIG. 2: Distribution of the counts obtained in one second on 
the vectors \VH),\HV),\XD),\DX),\RR), and \LL). Bars 
present the histograms of the measured numbers, and lines 
show the Poisson distribution with the mean values estimated 
from the experiment. Measurement time was 40 seconds for 
each vectors. 



where 



ty/2 - 2F 



merit time on a vector must be an integral multiple of 
one second. Figure [21 shows the histogram of the counts 
in one second on the vector 

B = {\VH), \HV), \XD), \DX), \RR), \LL)}, 

when the visibility of the two-photon states was esti- 
mated to be 0.92. The measurement time was 40 seconds 
on each vector. The distribution of the coincidence events 
obeys the Poisson distribution. Only small numbers of 
counts were observed on the vectors \HV) and \VH). 
Those observations agree with the prediction, therefore, 
we expect that the hypothesis testing in the previous sec- 
tions can be applied. 

In the following, we compare four testing meth- 
ods on experimental data with the fixed total time t. 
The testing method employ the different time allocations 

{tHH, tw, tDD, txX,tRL, t^R, tjIV,tVHi tnX,txD, tRR, thh} 

between the measurement vectors: 

(i) Modified visibility method: A is unknown. The 

coincidence and the anti-coincidence are measured 
with the equal time allocation; 

tHH = tw = tuD = txx = tRL = tLR 

t 

— tnv — tvH — tax — txD — tRR — tLL — t^. (26) 

(ii) Design I: A is unknown. The counts on coincidence 

and anti-coincidence are measured with the opti- 
mal time allocation at the target threshold Fq < 
0.899519; 



V2F + 1 + V2 - 2F 



t-2 



ty/2F + 1 



y/2F + 1 + V2 - 2F 



(28) 



(iii) Design II: A is known. Only the counts on anti- 
coincidence are measured with the equal time allo- 
cation at the target threshold Fo > 1/4; 

tHV = Wh = tDX = txD = tRR = tLL = 77 • (29) 

o 

(iv) Design III: A is known. Only the counts on anti- 
coincidence are measured. The time allocation is 
given by the two-stage method: 

tf 

tHV = tvH = tDX = txD = tRR, = tLL = ~T (30) 

6 

in the first stage, and 
m 



txy 



xy 



(t - tf) 



(31) 



tHH — tw — t'DD — txx 



tRL — tLR. — ~T 

6 



tHV — tvH — tDX — txD — tRR — tLL — ~Z , (27) 

6 



in the second stage. The observed count m xy in 
the first stage determines the time allocation in the 
second stage. 

We have compared the p-values at the fixed thresh- 
old Fo = 7/8 = 0.875 with the total measurement time 
t — 240 seconds. As shown in section III Bl the p-value 
measures the minimum risk probability to reject the hy- 
pothesis Hq, i.e., the probability to make an erroneous 
decision to accept insufficiently entangled states with the 
fidelity less than the threshold. The results of the exper- 
iment and the analysis of obtained data are described in 
the following. 
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In the method (i), we measured the counts on each 
vectors for 20 seconds. We obtained n\ — 9686 and n 2 — 
868 in the experiment, which yielded the p-value 0.343. 

In the method (ii), the optimal time allocation was 
calculated with l|28l) to be t\ — 55.6 seconds and t 2 — 
184.4 seconds. However, since the time allocation should 
be the integral multiple of second in our experiment, we 
used the time allocation t\ = 54 and t 2 = 186. That 
is, we measure the count on each coincidence vectors for 
9 seconds and on each anti-coincidence vectors for 31 
seconds. We obtained n\ = 7239 and n 2 = 2188 in the 
experiment, which yielded the p-value 0.0715. 

In the method (iii), we measured the count on each 
anti-coincidence vectors for 40 seconds. We used A = 
290 estimated from another experiment. We obtained 
n = 2808 in the experiment, which yielded the p-value 
0.0438. 

In the method (iv), the calculation is rather compli- 
cated. Similarly to (iii), A was estimated to be 290 from 
another experiment. In the first stage, we measured the 
count on each anti-coincidence vectors for tf/6 = 1 sec- 
ond. We obtained the counts 6,3,13,20,11, and 23 on 
the vectors \HV), \VH), \DX), \XD), \RR), and \LL), 
respectively. We made the time allocation of remain- 
ing 234 seconds for the second stage according to IpHjl. 
and obtained t HV = 28.14, t VH = 19.90, t DX = 41.42, 
txD = 51.37, tjm = 38.10, and tLL — 55.09. Since 
the time allocation should be the integral multiple of 
second in our experiment, we used the time allocation 
{t H v,tvH,tDX,txD,t RR ,t LL } = {28,20,42,51,38,55}. 
We obtained the counts on anti-coincidence nuv = 
99, uvh = 66, njjx = 703, uxd = 863, Urr = 531, and 
till = 853. Applying the counts and the time alloca- 
tion to the formula (|24|l . we obtained the p- values to be 
0.0310. 

The p-values obtained in the four methods are sum- 
marized in the table. We also calculated the p-values 
at different values of the threshold Fq as shown in Fig. 
13 We fixed time allocation for design I at t\ = 54 s 
and t 2 = 186 s. As clearly seen, the optimal time al- 
location between the coincidence vectors measurement 
and the anti-coincidence vectors measurement reduces 
the risk of a wrong decision on the fidelity (the p-value) 
in analyzing the experimental data. The counts on the 
anti-coincidence vectors is much more sensitive to the 
degradation of the entanglement. This matches our in- 
tuition that the deviation from zero provides a more ef- 
ficient measure than that from the maximum does. The 
comparison between (iii) and (iv) shows that the risk 
can be reduced further by the time allocation between 
the anti-coincidence vectors, as shown in Fig. [3] The 
optimal (Neyman) allocation implies that the measure- 
ment time should be allocated preferably to the vectors 
that yield more counts. Under the present experimental 
conditions, the optimal allocation reduces the risk proba- 
bility to about 75 %. The improvement should increased 
as the fidelity. However, the experiment showed almost 
no gain when the visibility was larger than 0.95. In such 



high visibility, errors from the maximally entangled state 
are covered by dark counts, which are independent of the 
setting of the measurement apparatus. 





(i) 


(ii) 


(iii) 


(iv) 


p-value at 0.875 


0.343 


0.0715 


0.0438 


0.0310 




0.872 0.874 0.876 0.878 

Threshold F 



FIG. 3: Calculated p-value as a function of the threshold. 
Dash-dot: (i)the modified visibility, dash: (ii)design I, dot: 
(iii)design II, solid: (iv) design III. 



IX. CONCLUSION 

We have applied the formulation of the hypothesis test- 
ing scheme and the design of experiment for the hypoth- 
esis testing of entanglement to the two-photon state gen- 
erated by SPDC. Using this scheme, we have handled 
the fluctuation in the experimental data properly. It has 
been experimentally demonstrated that the optimal time 
allocation improves the test in the terms of p-values: the 
measurement time should be allocated preferably to the 
anti-coincidence vectors in order to reduce the minimum 
risk probability. This design is particularly useful for 
the experimental test, because the optimal time alloca- 
tion depends only on the threshold of the test. We don't 
need any further information of the probability distri- 
bution and the tested state. We have also experimen- 
tally demonstrated that the test can be further improved 
by optimizing time allocation among the anti-coincidence 
vectors by using the two-stage method, when the error 
from the maximally entangled state is anisotropic. 
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